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Newton’s Law of Heating Simulation 
(Part 1: Constant Ambient Temperature) 


1 Introduction 


This simulation concerns itself with the temperature of a glass of iced tea, 
initially at 1°C, after it is placed in a refrigerator with a constant interior 
ambient temperature of 6°C. For this simulation, the object’s temperature is 
assumed to obey Newton’s Law of Heating. 


2 Mathematical modeling 
Starting with Newton’s Law of Heating, 
dT/dt = 1/2/min (6°C — T) with T(0) = 1°C (1) 
(Elected to use simple fractions instead of dealing with decimals.) 
Replacing dT /dt with T and omitting units for clarity, 
T = 1/2 (6-T) with T(0) =1 (2) 
Using differential equation techniques, 
T = -5e(-t/2) +6 (3) 


This analytical solution will be compared with the analog computation 
solution on page 3. 


3 Computer setup 
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Figure 1: Computer setup for Newton’s Law of Heating Simulation (Ta = constant) 


4 Results 





Figure 2: Temperature vs time simulation (constant ambient temperature)* 

















Analog Computation Analog Computation Analytical 

t(s) T(°C) T(°C) 
Oscilloscope Oscilloscope Estimate Calculator 
Scale: 0.5 min per small Scale: 0.4°C per small division Rounded to two 
division decimal places 
0.00 1.0 1.0 

5.00 5.6 5.6 

10.0 6.0 6.0 = Tmax 














Table 1: Solution Comparison: Analog Computation vs Analytical 


Newton’s Law of Heating Simulation 
(Part 2: Linearly Decreasing Ambient Temperature) 


1 Introduction 


This part of the simulation concerns itself with the temperature of a glass of 
iced tea, initially at 1°C, after it is placed in a refrigerator with an interior 
ambient temperature initially at 6°C but then begins to decrease linearly. For 


this simulation, the object’s temperature is assumed to obey Newton’s Law of 
Heating. 


2 Mathematical modeling 
Starting with Newton’s Law of Heating and Ta = (-1/3 °C/min t + 6°C), 
dT/dt = 1/2/min (-1/3°C/min t + 6°C — T) with T(0) = 1°C (1) 
(Elected to use simple fractions instead of dealing with decimals.) 
Replacing dT /dt with T and omitting units for clarity, 
T = 1/2 (-t/3 + 6-T) with T(0) =1 (2) 
Rewriting (2), 


T+T/2=-t/6+3 (3) 


First, seek a complementary solution (T-) by setting the right-hand side of (3) 
equal to zero. 


Tce +Tc/2 =0 
ToS te 2 
Integrating by inspection, 
In(T-) = -t/2 + constant 
Setting the constant = In(q), 
In(T./a) = -t/2 
T.= ae*(-t/2) (4) 
Next, seek a particular solution (T)) by selecting a trial function. 
Assuming Tp = Bt + y 
Tp=B (S) 
Substituting the expressions for (4) and (5) into (3), 
B+(Bt+ y)/2 =-t/6+ 3 
B/2t+B+y/2=-t/6+3 (6) 
Comparing coefficients of like terms on each side of (6), 
B/2 =-1/6 
B=-1/3 
and 


B+y/2=3 


-1/3 + y/2 =3 
y=20/3 
Thus, 
Tp = -t/3 + 20/3 
and 
T=Tce+Tp=ae(-t/2) —t/3 + 20/3 
Invoking the initial condition T(0) = 1, 
1=a+20/3 
a= -17/3 
yielding 
T = -17/3 e(-t/2) —t/3 + 20/3 (6) 


This analytical solution will be compared with the analog computation 
solution on page 7. 


Differentiating (3) and setting it equal to zero will yield Tmax. 
T = 17/6 e*(-t/2) — 1/3 
0=17/6e%(-t’/2)-1/3 

Using algebra, 

t’ = -2In(2/17) = 4.28 
and 


Tmax & -17/3 e*(-4.28/2) — 4.28/3 + 20/3 = 4.57 


3 Ambient temperature generator computer setup 
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Ta = -t/3 + 6 
1 Th=t/3-6 


Figure 1: Computer setup to generate decreasing ambient temperature 


4 Final computer setup 
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Figure 2: Computer setup for Newton’s Law of Heating Simulation 


5 Results 
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Figure 3: Temperature vs time simulation (decreasing ambient temperature)* 























Analog Computation Analog Computation Analytical 
t(s) T(°C) T(°C) 
Oscilloscope Oscilloscope Estimate Calculator 
Scale: 0.5 min per small division | Scale: 0.4°C per small division 

0.00 1.0 1.0 

10 3.6 3.3 





Table 1: Solution Comparison: Analog Computation vs Analytical 


*For this application note, each display was produced during a single run by a 
differential equation analog computer prototype using discrete components 
with tolerances between 1% and 10%. 
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